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(Dated: March 30, 2007) 

We present a study of events with Z bosons and associated jets produced at the Fermilab Tevatron 
Collider in pp collisions at a center of mass energy of 1.96 TeV. The data sample consists of nearly 

14,000 Z/^* e^e^ candidates corresponding to an integrated luminosity of 0.4 fb~^ collected 
with the D0 detector. Ratios of the Z/^*+ > n jet cross sections to the total inclusive Z/^* cross 
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section have been measured for n = 1 to 4 jets, and found to be in good agreement with a next-to- 
leading order QCD calculation and with a tree-level QCD prediction with parton shower simulation 
and hadronization. 

PACS numbers: 13.38.Dg, 14.70.Hp, 13.87.-a, 12. 38. Aw, 12.38.Qk, 13.85.-t 



Leptonic decays of electroweak gauge bosons, W and 
Z, produced in association with jets are prominent sig- 
natures at present and future hadron colliders. Measure- 
ments of W (or Z) + > n jet cross sections are important 
for understanding perturbative quantum chromodynam- 
ics (QCD) calculations and for developing Monte Carlo 
(MC) simulation programs capable of handling partons 
in the final state at leading order (LO), or, in some cases, 
next-to- leading order (NLO). Furthermore, the produc- 
tion of or Z bosons with associated jets represents a 
significant background to Higgs boson searches, as well as 
to other standard model processes of interest, such as top 
quark production, and many searches for new phenom- 
ena at the Fermilab Tevatron Collider and at the CERN 
Large Hadron Collider. 

Measurements of Z -I- > n jet cross sections with lower 
integrated luminosity and at lower center of mass energy 
were performed previously by the CDF collaboration [l| . 
In this study, we present the first measurement of the 
fully corrected ratios of the Z/'^*+ > n jet production 
cross sections to the total inclusive Z/j* cross section 
for jet multiplicities n — 1 — 4 in collisions at ^/s — 
1.96 TeV. Cross section measurements based on inclusive 
jet multiplicities provide theoretically sound observables, 
and can be compared to a variety of predictions. Our 
results are based on a data sample corresponding to an 
integrated luminosity of 0.4 fb^^ accumulated with the 
D0 detector. 

The elements of the D0 detector [2[ of primary im- 
portance to this analysis are the uranium/liquid-argon 
sampling calorimeter and the tracking system. The D0 
calorimeter has a granularity of At; x Acj) = 0.1 x 0.1, 
forming projective towers, where rj is the pseudorapid- 
ity (^ = — ln[tan(0/2)], 9 is the polar angle relative to 
the proton beam), and (j) is the azimuthal angle. The 
calorimeter has a central section covering pseudorapidi- 
ties up to « 1.1, and two end calorimeters that extend 
the coverage to |?7| ~ 4.2. The tracking system consists 
of a silicon micro-strip tracker and a central fiber tracker, 
both located within a 2 T superconducting solenoidal 
magnet, with designs optimized for tracking and vertex- 
ing at pseudorapidities of \ri\ < 3 and \r]\ < 2.5, respec- 
tively. 

The data sample for this analysis [3] was collected be- 
tween April 2002 and June 2004. Events from Z/7* 
e+e" decays were selected with a combination of single- 
electron triggers, based on energy deposited in calorime- 
ter towers (Arj x Acj) = 0.2 x 0.2). Final event selection 
was based on detector performance, event properties, and 
electron and jet identification criteria. 



Events were required to have a reconstructed pri- 
mary vertex with a position along the beam direction 
within 60 cm of the detector center. Electrons were re- 
constructed from electromagnetic (EM) clusters in the 
calorimeter using a simple cone algorithm. The two elec- 
tron candidates in the event with the highest transverse 
momentum components relative to the beam direction 
(pt), and both with pr > 25 GeV, were used to recon- 
struct the Z boson candidate. The two electrons were 
required to be in the central region of the calorimeter 
I'ydotI < 1-1 (pseudorapidity 77dct is calculated relative to 
the center of the detector), and at least one required to 
fire the trigger(s) for the event. The electron pair also 
had to have an invariant mass consistent with the Z bo- 
son mass of 75 GeV < Mee < 105 GeV. 

To reduce background (mainly from jets misidentified 
as electrons), the EM clusters were required to pass three 
quality criteria based on the shower profile: (i) the elec- 
tron had to deposit at least 90% of its energy in the 
2 1-radiation- length EM calorimeter (ii) the lateral and 
longitudinal shape of the energy cluster had to be consis- 
tent with those of an electron, and (iii) the electron had to 
be isolated from other energy deposits in the calorimeter, 
with an isolation fraction /iso < 0.15. (The isolation frac- 
tion is defined as /iso = [-^(0.4) - Sem(0.2)]/£;em(0.2), 
where -E(i?conc) and -E'EM(J^conc) are respectively the to- 
tal and EM energies within a cone of radius i?conc — 
\J (At])^ -I- (A0)^ centered around the direction of the 
electron.) Additionally, at least one of the electrons was 
required to have a spatially matched track associated 
with the reconstructed calorimeter cluster, and the track 
momentum had to be consistent with the energy of the 
EM cluster. A total of 13,893 events passed the selection 
criteria. 

Jets were reconstructed using the "Run II cone al- 
gorithm" Q that combines particles within a cone of 
radius i?conc — 0.5. Spurious jets from isolated noisy 
calorimeter cells were eliminated through selections on 
patterns of jet energy deposition. Jets were required 
to be consistent with energy depositions measured at 
the trigger stage. This requirement was introduced to 
address precision readout noise problems: The jet en- 
ergy at the Level 1 trigger tower level was compared 
to the jet energy derived from the jet cone algorithm, 
which was based on calorimeter cell precision readout. 
The transverse momentum of each jet was corrected for 
multiple pp interactions, calorimeter noise, out-of~cone 
showering effects, and energy response of the calorimeter 
as determined from the missing transverse energy bal- 
ance of photon-jet events Q. Jets were required to have 
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Pt > 20 GeV and \'q\ < 2.5, and were eliminated if they 
overlapped with electrons from Z boson decay within 
Ai? = ^(A77)2 + (A0)2 = 0.4. Small losses of jets re- 
sulting from this separation criterion for electrons from 
Z boson decays were estimated as a function of the num- 
ber of associated jets using a pythia Q MC sample. 

The jet energy resolutions were derived from a mea- 
surement in photon-l-jet data for low jet energies and dijet 
data for higher jet energy values. Fits to the transverse 
energy asymetry [pt(1) — Pt(2)]/[pt(1) +Pt(2)] between 
the transverse momenta of the back-to-back jets and/or 
photon {pt(1) and Pt(2)) were then used to obtain the 
jet energy resolution as a function of jet rapidity and 
transverse energy. The largest contribution to the jet en- 
ergy resolution uncertainty was due to limited statistics 
in the samples used. 

The electron efficiencies for trigger, track matching, 
reconstruction, and identification were determined from 
data, based on a "tag-and-probe" method. Z candi- 
dates were selected with one electron (the tag) satisfy- 
ing a tighter track-matching requirement to further re- 
duce background contamination, and another electron 
(the probe) with all other criteria applied, except the 
one under study. The fraction of events with probe elec- 
trons passing the requirement under study determined 
the efficiency of a given criterion. The overall trigger 
efficiency for Z candidates that survived the analysis se- 
lections was found to be greater than 99%. The electron 
reconstruction and identification efficiencies were mea- 
sured as a function of azimuthal angle and pT^ and the 
average efficiency was found to be about 89%. The com- 
bined spatial and energy track-matching efficiency was 
measured to be about 77%. The electron reconstruction, 
selection, trigger, and track-matching efficiencies were 
examined as a function of jet multiplicity. No signifi- 
cant variations of the efficiencies were observed, except 
for the track-matching efficiency, for which the multiplic- 
ity dependence was taken into account in correcting the 
data. 

The kinematic and geometric acceptance for electrons 
from Z/ 7* decays in the mass region of 75 GeV < Mee < 
105 GeV, for a primary vertex within 60 cm of the de- 
tector center, was determined as a function of jet multi- 
plicity. An inclusive pythia sample was used to calcu- 
late the acceptance for the inclusive Z/j* sample. The 
PYTHIA events were weighted so that the pr distribu- 
tion of the Z boson in the MC agreed with data. The 
jet multiplicity dependence of the acceptance was calcu- 
lated using a Z/"/* +n parton leading-order generator Q, 
with the evolution of partons into hadrons carried out in 
PYTHIA. All the samples were processed through full D0 
detector simulation using GEANT Q and the D0 recon- 
struction software. The overall dielectron acceptance for 
the Z/j*+ > 4 jet sample was found to be about 30% 
higher than the acceptance for the Z/"f* inclusive sample, 
because events with jets tend to recoil from Z bosons of 



larger px, and thereby produce decay products that are 
more likely to fall within the geometric acceptance. 

The reconstruction and identification efficiency of jets 
was determined from a MC sample with full detector sim- 
ulation, and processed through the same programs as the 
data. A scaling factor was applied to the MC jets to ad- 
just their reconstruction and identification efficiency to 
that of jets in data, using the "Z pr-balance" method [§]. 
In events with Z candidates, a search was performed for 
a recoiling jet opposite in azimuth to the Z boson. The 
probability of finding a recoiling jet as a function of the 
Pt of the Z was measured in data and MC. The ratio 
of these probabilities defined the scaling factor that was 
applied to the MC jets. After applying the scaling fac- 
tor, the jet reconstruction and identification efficiency 
was determined by matching particle-level jets (i.e., jets 
found from final-state generator-level particles, after par- 
ton hadronization) to calorimeter jets. The efficiency was 
parameterized as a function of the px of the particle-level 
jet, where the pr values were smeared with jet energy 
resolutions observed in data, as measured in three ij re- 
gions of the calorimeter. As a cross check, the scaling 
factor determined from the "Z p^-balance" method was 
compared to the scaling factor obtained for photon-|-jet 
events, and found to be consistent with one another. 

The primary background to the Z/j* dielectron sig- 
nal is from multijet production, in which the jets have 
a large electromagnetic component or they are mismea- 
sured in some way that causes them to pass the electron 
selection criteria. We refer to this instrumental back- 
ground as "QCD". For the Z/j*+ > — 2 jet samples, 
a convoluted Gaussian/Breit-Wigner function was used 
to fit the Z lineshape, and an exponential form was used 
to account for both the QCD background and the Drell- 
Yan (7*) component of the signal. For the lower statis- 
tics Z/j*+ > 3 jet sample, the contributions from QCD 
and Droll- Yan components were estimated from the side 
bands of the Z in the dielectron invariant mass spectrum. 
In each case, a pythia sample was used to disentangle 
the QCD component from the Drell-Yan contribution. 
The background contribution for the Z/'y*+ > 4 jet mul- 
tiplicity sample was estimated by extrapolating to n = 4 
an exponential fit to the QCD background in the — 3 
jet multiplicity bins. The background contribution from 
QCD processes was found to be 3 — 5%, depending on 
jet multiplicity. There are also contributions to Z/j* 
candidates that are not from misidentified electrons, but 
correspond to other standard model processes (e.g., ti 
production, Z t+t~, W ev). These small (< 1%) 
irreducible background contributions were also taken into 
account in the analysis. 

The cross sections as a function of jet multiplicity were 
corrected for jet reconstruction and identification effi- 
ciencies, and for event migration due to the finite jet 
energy resolution of the detector. The correction fac- 
tors were determined using two independent MC sam- 
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TABLE I: Cross-section ratios (Rn) with statistical and systematic uncertainties (all xlO ^) for different inclusive jet multi- 
plicities. 



Multiplicity {Z/j*+ > n jets) 


n > f 


n>2 


n > 3 


n > 4 


Rn 


f20.f 


18.6 


2.8 


0.90 


Total Statistical Uncertainty 


±3.3 


±f.4 


±0.56 


±0.44 


Total Systematic Uncertainty 




— o.U, -\-v.Z 


— l.UD, +1.4o 


n /in _i_n a s 


Jet Energy Calibration 


iff. 7 


±3.3 


±0.74 


±0.23 


Jet Reconstruction/Identification 


-7.0, +2.2 


-2.9, +4.3 


-0.64, +0.82 


-0.30, +0.40 


Unsmearing Procedure 


-3.6, +2.2 


-f.6,+2.4 


-0.24, +0.85 


-0.08, +0.09 


Jet Energy Resolution 


-2.7, +3.4 


-0.04, +0.f3 


-0.f7, +0.f5 


-0.03, +0.04 


Acceptance 


±f.8 


±0.7 


±0.fO 


±0.003 


Efficiencies (Trigger, EM, Track) 


±8.5 


±f.3 


±0.20 


±0.07 


Electron- Jet-Overlap 


±3.2 


±0.7 


±0.f4 


±0.05 



pies, both tuned to match the measured inclusive jet 
multiplicity and jet pT distributions in data. The first 
sample was based on PYTHIA simulations. The second 
sample (ME-PS) was based on madgraph [l^ Z/7* + n 
LO Matrix Element (ME) predictions, using pythia for 
parton showe ring (PS) and hadronization, and a mod- 
ified CKKW [III method to map the Z/j* + n parton 
event into a parton shower history [T2I. The ME-PS 
predictions relied on MADGRAPH tree-level processes of 
up to three partons. Both these samples contained only 
particle-level jets (i.e., no detector simulation). The pr 
of the jets was smeared with the jet energy resolutions 
found in data. Subsequently, some jets were removed 
randomly from the sample, to simulate the measured 
jet reconstruction/identification efficiencies. The ratio 
of the two inclusive jet multiplicity distributions (the 
generated distribution and the one with the jet recon- 
struction/identification efficiency and energy resolution 
applied) determined the unsmearing correction factors 
for a given MC sample. The weighted averages of the 
correction factors corresponding to the two sets of MC 
procedures were applied as a function of jet multiplic- 
ity to correct the jet multiplicity spectrum in data. The 
differences between the correction factors for the two cal- 
culations contribute to the systematic uncertainty of the 
procedure. Another source of systematic uncertainty was 
determined from a closure test, and was estimated by 
applying the full unsmearing procedure to a MC control 
sample. The unsmearing correction factors range from 
f .11 to 2.9 for > 1 and > 4 jets, respectively. 

The fully corrected ratios, i?„, of the Z/j*+ > n jet 
production cross sections to the inclusive Z/ 7* cross sec- 
tion 



ct(Z/7*+ > n jets) 

<y{z/r) 



(1) 



for the mass region 75 GeV < Mee < 105 GeV are sum- 
marized in Table m Systematic uncertainties include con- 
tributions from jet energy calibration corrections, jet re- 
construction and identification efficiency, the unsmear- 
ing procedure, jet energy resolution, and variations in 
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FIG. 1: Ratios of the Z/'y*-\- > n jet cross sections to the 
total inclusive Z/j* cross section versus jet multiplicity. The 
uncertainties on the data (dark circles) include the combined 
statistical and systematic uncertainties added in quadrature. 
The dashed line represents predictions of LO Matrix Element 
(ME) calculations using pythia for parton showering (PS) 
and hadronization, normalized to the measured Z/'y*+ > 1 
jet cross-section ratio. The dotted line represents the predic- 
tions of PYTHIA normalized to the measured Z/^*+ > 1 jet 
cross-section ratio. The two open diamonds represent predic- 
tions from MCFM. 



the acceptance for different event generators. They also 
take into account uncertainties in the variation of effi- 
ciencies for the trigger, electron reconstruction, identifi- 
cation, and track matching as a function of jet multiplic- 
ity, as well as uncertainties due to the electron-jet overlap 
correction. All these uncertainties are assumed to be un- 
correlated, and are added in quadrature to estimate the 
total systematic uncertainty. The statistical uncertain- 
ties include contributions from the number of candidate 
events, background estimation, acceptance, efficiencies, 
and the unsmearing correction. 

Figure [T] shows the fully corrected measured cross- 
section ratios for Z/j*+ > n jets as a function of 
jet multiplicity, compared to three QCD predictions. 
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FIG. 2; Comparison between data and theory (ME-PS) for 
the highest pr jet distribution in the Z/y*+ > 1 jet sample 
(dark circles), for the second highest pr jet distribution in 
the Z/-y*+ > 2 jet sample (open circles), and for the third 
highest Pt jet distribution in the Z/'y*+ > 3 jet sample (open 
triangles). The uncertainties on the data are only statistical. 
The MC distributions are normalized to the data. 



ratios were found to be in good agreement with MCFM 
and an enhanced leading-order matrix element prediction 
with PYTHiA-simulated parton showering and hadroniza- 
tion. PYTHIA simulations alone appear to exhibit a deficit 
in high jet multiplicity events. 

We thank S. Mrenna for providing us the ME-PS MC 
sample. We thank the staffs at Fermilab and collaborat- 
ing institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); FASI, 
Rosatom and RFBR (Russia); CAPES, CNPq, PAPER J, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); POM (The Netherlands); PPARC (United 
Kingdom); MSMT (Czech Repubhc); CRC Program, 
CFI, NSERC and WestGrid Project (Canada); BMBF 
and DFG (Germany); SFI (Ireland); The Swedish Re- 
search Council (Sweden); Research Corporation; Alexan- 
der von Humboldt Foundation; and the Marie Curie Pro- 
gram. 



MCFM is a NLO calculation for up to Z/j* + 2 par- 
ton processes. CTEQ6M lj| parton distribution func- 
tions (PDF) were used in MCFM, and the factorization 
and rcnormalization scales fxp, fJ-R were both set to the 
Z boson mass, Mz- Varying the PDF set and the renor- 
malization/factorization scales to +Px z had a mini- 
mal effect on the mcfm cross-section ratios. The ME-PS 
predictions are normalized to the measured Z/"f*+ > 1 
jet cross-section ratio, and use the CTEQ6L PDF, with 
the factorization scale set to fj^p = Mz, and the rcnor- 
malization scale set to ^.r = pTjet for jets from initial 
state radiation and fiji = kxjet for jets from final state 
radiation {kxjet is the transverse momentum of a radi- 
ated jet relative to its parent parton momentum). The 
PYTHIA predictions are also normalized to the measured 
^/7*+ > 1 jet cross-section ratio. Here, CTEQ5L [15 1 
PDFs are used, and the factorization and rcnormaliza- 
tion scales are set to — f^R = Mz- The mcfm and 
ME-PS predictions are generally in good agreement with 
the data, pythia predicts fewer events at high jet mul- 
tiplicity because of missing higher order contributions at 
the hard-scatter level. 

Figure [5] compares jet pT spectra of the n*'' jet, n = 
1,2,3, in Z/7*-|- > n jet events to the ME-PS MC predic- 
tions. The MC events have been passed through the full 
detector simulation, and the jet pr spectra normalized 
separately to the data distributions. Good agreement 
can be seen over a wide range of jet transverse momenta. 

In summary, we have presented the first measurements 
of fully corrected ratios of the Z/^*+ > n jet (n = 1 — 4) 
production cross sections to the total inclusive Z 1^* cross 
section in pp collisions at ^/s — 1.96 TeV. The measured 
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